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Abstract. We present extensive evolutionary models of stars 
with initial zero-metallicity, covering a large range of initial 
masses (i.e. 0.7M© < M < 1OOM ). Calculations are carried 
out at constant mass, with updated input physics, and applying 
an overshooting scheme to convective boundaries. The nuclear 
network includes all the important reactions of the p-p chain, 
CNO-cycle and a-captures, and is solved by means of a suit- 
able semi-implicit method. The evolution is followed up to the 
thermally pulsing AGB in the case of low- and intermediate- 
mass stars, or to the onset of carbon burning in massive stars. 

The main evolutionary features of these models are dis- 
cussed, also in comparison with models of non-zero metallic- 
ity. Among several interesting aspects, particular attention has 
been paid to describe: i) the first synthesis of 12 C inside the 
stars, that may suddenly trigger the CNO-cycle causing partic- 
ular evolutionary features; ii) the pollution of the stellar sur- 
face by the dredge-up events, that are effective only within par- 
ticular mass ranges; iii) the mass limits which conventionally 
define the classes of low-, intermediate-, and high-mass stars 
on the basis of common evolutionary properties, including the 
upper mass limit for the achievement of super-Eddington lumi- 
nosities before C-ignition in the high-mass regime; and iv) the 
expected pulsational properties of zero-metallicity stars. 

All relevant information referring to the evolutionary tracks 
and isochrones is made available in computer-readable format. 

Key words: stars: evolution - stars: interiors - stars: Hertz- 
sprung-Russell (HR) diagram - stars: low-mass 



1. Introduction 

The standard Big Bang Nucleosynthesis (e.g. Olive 1999) pre- 
dicts that essentially no elements heavier than 7 Li should exist 
in the gas mixture in which the first stars are generated. Only 
in the framework of inhomogeneous Big Bang nucleosynthe- 
sis (see Jedamzik 2000 and references therein) we can expect 
some non-negligible primordial metal production. 

The stars with negligible, if not zero, initial metal abun- 
dance - named population III objects (Pop-Ill for short) - are 
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expected to have structure and evolution distinct from those in 
which traces of heavier elements, like those of the CNO group, 
are present. 

The major difference between zero-metal stars and those 
of normal metal content (even down to very low values of Z) 
lies in the mechanism of nuclear energy generation. In fact, 
owing to the lack of CNO nuclei, the pre-main sequence grav- 
itational contraction cannot stop until the central temperature 
and density are high enough to allow the p-p chain to provide 
the energy budget. Since the p-p chain is a poor thermostat as 
compared to the CNO-cycle, very high temperatures can be 
reached in the central regions. This may eventually lead (de- 
pending on the mass of the star) to the first synthesis of pri- 
mary carbon through the 3-a reaction, while the star is still 
on the main sequence. As a consequence, the CNO-cycle is 
activated at the very high temperatures characterising the 3-a 
process, possibly causing a dramatic change in the dominant 
energy source (from p-p chain to CNO-cycle), hence affect- 
ing the stellar structure amid the main sequence evolution. As 
demonstrated below, similar changes in the evolutionary be- 
haviour may also occur at later stages, whenever carbon is first 
produced (or transported) into different regions of the star. 

The threshold abundance of 12 C at which the distinct be- 
haviour of zero metal stars appears is as low as Z w Ac ~ 
10~ 10 - 10" 9 (see Cassisi & Castellani 1993). Starting from 
these values and higher, the efficiency of CNO-cycle is suffi- 
cient to recover the standard behaviour. 

Needless to say, the distinct evolutionary behaviour of these 
stars may also imply a very distinct nucleosynthesis and chem- 
ical pollution of the interstellar medium, if compared to stars 
with initial Z > 10~ 10 . Also their properties in the HR- 
diagram should be quite peculiar. Both aspects are particularly 
important in the modelling of 

1 . the initial chemical enrichment of galaxies, when the first 
metals and additional light elements have been made avail- 
able to the second generation of stars, and probably also to 
the intergalactic medium; 

2. the initial spectrophotometric evolution of galaxies, when 
the first ionizing photons of stellar origin probably con- 
tributed to the reionization of the Universe; 

3. the observational characteristics that the remaining low- 
mass Pop-Ill stars (if any were formed) would present to- 
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day, in order to guide observational searches for these ob- 
jects. 

In this context, the present paper is aimed to provide the 
necessary tools for the investigation of the above-mentioned 
questions. We start by presenting an extensive and homoge- 
neous set of evolutionary tracks for Pop-Ill stars, covering the 
complete mass range from 0.7 to 100 Mq. All of the models 
are computed with convective overshooting and under the as- 
sumption of constant-mass evolution. A subsequent paper will 
deal with the case of more massive stars (M > 100 Mq), in 
which stellar winds had likely been driven by a state of pulsa- 
tion instability. 

The plan of the paper is as follows. After briefly reviewing 
the available literature on evolutionary models of primordial 
stars (Sect. ||), the main assumptions of present calculations are 
described in Sect. [|concerning the nucleosynthesis, and Sect. ^ 
for the other physics prescriptions. Then, the presentation and 
discussion of the results start with Sect. [| - dealing with gen- 
eral evolutionary aspects (e.g. the onset of the 3-a reaction and 
the activation of the CNO cycle, lifetimes, etc.) - and proceed 
through Sects. || and [|, which focus on particular issues of 
low- and intermediate-/high- mass models, respectively (e.g. 
the AGB phase, the Eddington critical luminosity). Our predic- 
tions for possible changes in the surface chemical composition 
are given in Sect. ||. Finally, Sect. ^ is dedicated to an analysis 
of the expected pulsational properties of Z = stars with dif- 
ferent masses and evolutionary stages, to investigate whether 
pulsation might be a discriminating feature in observational 
searches of Pop-Ill objects. Appendices describe the tables of 
the Z = evolutionary tracks and related isochrones available 
in electronic format. 

2. A quick review of past work on Pop-Ill stars 

The mass-scale of Pop-Hi stars is still a matter of debate. The 
once widely accepted idea that only massive, very massive or 
even super-massive stars could form in the primordial medium 
has long limited the study of Pop-Ill stars to objects more mas- 
sive than say 10 Mq, possibly in the range 10 2 to 10 5 Mq 
(Carr et al. 1982). Similar suggestions were made by Hutchins 
(1976), Silk (1977) and Tohline (1980) whose work indicated 
a minimum mass of 150 — 200 Mq. 

However, over the years, interest has turned towards more 
"normal" stars. In fact, various fragmentation models of pri- 
mordial gas have suggested that the first stars might have 
formed with low/intermediate masses, e.g. Yoshi & Saio (1986) 
find the peak of the mass function at roughly 4—10 M Q ; Naka- 
mura & Umemura (1999) find a the typical mass of 3M , 
which may grow to 16 Mq by accretion; and Nakamura & 
Umemura (2000a) find a lower mass-limit extending down to 
1 Mq. With more detailed and extended calculations for the 
collapse and fragmentation of primordial gas clouds, Naka- 
mura & Umemura (2000b) suggest that the initial mass func- 
tion of Pop-Ill stars may be bimodal with peaks at about 1 Mq 
and 100 Mq, the relative heights of these peaks being a func- 
tion of the collapse epoch. 



The provisional conclusion we may get from all this is that 
the mass peak of Pop-Ill stars is likely greater than say 1 Mq 
and possibly in the range 1 to 10 Mq. However, both lower and 
much higher values cannot be firmly excluded. 

Ezer & Cameron (1971) first studied models of Pop-Ill 
stars in the mass range 5 to 200 Mq from the pre-main se- 
quence until core H-exhaustion. In a subsequent paper Ezer 
(1972) investigated the evolution of a 3 Mq star up to the stage 
of core He-exhaustion. Cary (1974) and Castellani & Paolicchi 
(1975) computed zero-age-main-sequence (ZAMS) models in 
the mass range 2 to 20 Mq and 1 to 8 Mq, respectively. 

Woosley & Weaver (1981) followed the evolution at con- 
stant mass of a 200 Mq star until the supernova explosion 
stage, thus allowing the first estimate of the contamination 
power of this type of star. Very massive Pop-Ill stars were first 
investigated by Bond et al. (1982, 1983). El Eid et al. (1983) 
and Ober et al. (1983) explored the mass range from 80 to 
500 Mq, including mass loss by stellar winds up to the stage of 
supernova explosion or collapse to a black hole. Similar studies 
were undertaken by Klapp (1983, 1984) for very massive stars 
in the range 500 to 10 4 Mq. These papers made it possible to 
get an estimate of the chemical yields both in the wind and ex- 
plosive stages. Forieri (1982) studied the evolution of massive 
Pop-Ill stars in mass range 10 to 100 Mq up to the central He- 
exhaustion stage. These models were evolved at constant mass 
but included convective overshoot. Castellani et al. (1983) pre- 
sented models for 10, 15 and 25 Mq stars from ZAMS to C- 
ignition. 

The advanced phases of low-mass stars remained almost 
unexplored till the early study by D'Antona (1982) of a 1 Mq 
star up to the stage of core He-flash, where for the first time the 
much less extended RGB was noticed (see also Castellani & 
Paolicchi 1975 for an earlier suggestion). Chieffi & Tornambe 
(1984) investigated the evolution of a 5 Mq object, noticing the 
absence of thermal instabilities during the double-shell (AGB) 
phase. Fujimoto et al. (1984) pointed out that threshold values 
for the core mass and CNO abundances exist for the occurrence 
of He-shell flashes during the AGB phase. Cassisi & Castel- 
lani (1993) presented grids of models with masses in the range 
0.7 < M < 15M and metallicity from Z = 10~ 4 down to 
Z = 10~ 10 . The models were followed from ZAMS till He- 
flash, final cooling of the CO core or C-ignition as appropriate. 
Cassisi et al. (1996) analysed the RR Lyrae properties and the 
unusually strong He-shell flashes of low-mass stars with very 
low initial metal content (—10 < log Z < —5). 

Many recent/ongoing works on zero-metallicity stars have 
been recently presented in the symposium proceedings edited 
by Weiss et al. (2000a), to which the reader should refer for 
a comprehensive picture of the present state of understanding 
of these objects. Particularly impressive are the studies being 
carried out by Heger et al. (2000) and Umeda et al. (2000), 
where massive zero-metallicity stars are made to evolve up to 
the pre-supernova stage and supernova explosion, respectively. 
These studies describe in detail the nucleosynthesis of heavy 
elements occurring in these stars. 
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In the work by Weiss et al. (2000b) the authors investigate 
the role of diffusion and external pollution in the evolution of 
metal-free low-mass stars, finding that no sizable effects oc- 
cur (the accreted metals cannot reach the central H-burning 
regions, which would have otherwise affected the structure 
of zero-metallicity stars). Using an extended nuclear network, 
they also demonstrate that hot p-p chains are less effective than 
the 3-a process in producing carbon inside a zero-metallicity 
star. 

Finally, we should mention the work by Schlattl et al. 
(2001), who investigate the possible changes in the envelope 
chemical composition as the result of convective dredge-up 
during the He-flash of Pop-Ill RGB stars (see also Fujimoto 
et al. 2000). 



3. Nuclear network and nucleosynthesis 

In general, the evaluation of temporal changes in the chemical 
abundances due to nuclear reactions requires the integration of 
the system of equations below: 

jY t = -[ij] + H Y r Y sl i = 1, ...N el (1) 

where N e \ is the number of elemental species involved; Y. L = 
Xi/Ai, Xi, and Ai denote the abundance by number (mole 
g _1 ), the abundance by mass, and the atomic mass of the el- 
emental species i, respectively. In the right-hand side of the 
equation (|l]), [ij] stands for the rate of the generic reaction 
which converts the element i into another nucleus because of 
the interaction with the element j, whereas [rs] is the rate of the 
generic reaction transforming nuclei s and r into the element i. 

In this work the adopted nuclear network consists of the 
sets of reactions for the ppl -chain 
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and the most important a-capture reactions 
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The abundance equations (jl|) are solved with the aid of a 
semi-implicit extrapolation scheme (Bader & Deuflhard 1983), 
which combines all the considered H- and He-burning reac- 
tions (Eqs. ||, ||, and without any assumption for nuclear 
equilibria. Such procedure is necessary to follow the nucle- 
osynthesis occurring in stars with original zero (or extremely 
low) metal content, for various reasons. 

The simultaneous consideration of H- and He-burning re- 
actions is required by the fact that in these stars H-burning may 
occur at temperatures and densities high enough to allow the 
synthesis, via a-captures, of the missing catalysts (i.e. 12 C and 
le O) for the activation of the CNO cycle. 

The choice of a semi-implicit method is motivated by the 
search for a convenient compromise between the higher accu- 
racy typical of the explicit scheme, and the better stability of 
the solution guaranteed by the implicit scheme. 

We recall that, from a technical point of view, in im- 
plicit/explicit methods the increments of the dependent vari- 
ables over a given integration step are calculated from the 
derivatives evaluated at the new/old location of the indepen- 
dent variable. The net effect is that implicit methods are stable 
for large integration steps, but with a certain loss of accuracy 
in following the solution towards the equilibrium. The reverse 
situation occurs for explicit methods. 

Actually, implicit methods are required when dealing with 
stiff differential equations, in which the changes of the depen- 
dent variables involve quite different scales of the independent 
variable. This feature matches the case of the abundance equa- 
tions, as different elemental species are characterised by very 
different nuclear lifetimes. For instance, at temperatures near 
25 x 10 6 K, the e-folding time for 15 N is on the order of years, 
whereas for 14 N it is on the order of 10 5 years. In general, im- 
plicit methods converge to the true equilibrium solutions of the 
nuclear network for relatively large time steps, whereas purely 
explicit methods would require extremely short time steps to 
both find the solution and maintain stability. Alternatively, in 
order to make the explicit scheme less time consuming this 
may be coupled with external assumptions on nuclear equilib- 
ria (usually made for the CNO nuclei). This possible choice 
is however quite risky, specially in the case of zero-metallicity 
stars in which the total abundance by number of the CNO el- 
ements during H-burning is not constant, but rather increasing 
thanks to the contribution from He-burning. 



4 



P. Marigo et al.: Zero-metallicity stars 



4. General physical input of the stellar models 

Chemical composition. The adopted initial composition con- 
sists of a mixture of just hydrogen and helium, with mass frac- 
tions of X = 0.77 and Y = 0.23, respectively. The adopted 
Y value is consistent with most recent determinations of the 
primordial helium content, that typically range from 0.230 to 
0.234 (see e.g. Peimbert 1996 for a review; or Olive et al. 1997 
and Viegas et al. 2000 for recent determinations). 

Opacity. The radiative opacities are from the OPAL group 
(Rogers & Iglesias 1992; Iglesias & Rogers 1993) for tem- 
peratures higher than logT = 4.1, and from Alexander & 
Ferguson (1994) for logT < 4.0. In the temperature interval 
4.0 < logT < 4.1, a linear interpolation between the opacities 
derived from both sources is adopted. The conductive opaci- 
ties of electron-degenerate matter are from Hubbard & Lampe 
(1969). 

Equation of state. The equation of state (EOS) for tempera- 
tures higher than 10 7 K is that of a fully-ionized gas, including 
electron degeneracy in the way described by Kippenhahn et 
al. (1965). The effect of Coulomb interactions between the gas 
particles at high densities is introduced following the prescrip- 
tion by Straniero (1988; cf. also Girardi et al. 1996). 

Nuclear reaction rates and neutrino losses. The reaction rates 
are from the compilation of Caughlan & Fowler (1988), ex- 
cept for 17 0(p, a) 14 N and 17 0(p, 7) 18 F, for which we use the 
more recent determinations by Landre et al. (1990). The un- 
certain 12 C(a, 7) le O rate was set to 1.7 times the values given 
by Caughlan & Fowler (1988), as indicated by the study of 
Weaver & Woosley (1993) on the nucleosynthesis by massive 
stars. The electron screening factors for all reactions are those 
from Graboske et al. (1973). 

The energy losses by pair, plasma, and bremsstrahlung neu- 
trinos, are from Haft et al. (1994). 

External convection and overshooting. The energy transport in 
the outer convection zone is described according to the mixing- 
length theory of Bohm-Vitense (1958). The mixing length pa- 
rameter a = 1.68 comes from the calibration of a solar model 
(Girardi et al. 2000). 

The extension of convective boundaries is estimated by 
means of an algorithm which takes into account overshooting 
from both core and envelope convective zones. The formal- 
ism is fully described in Bressan et al. (1981) and Alongi et 
al. (1991). The main parameter describing overshooting is its 
extent A c across the border of the convective zone, expressed 
in units of pressure scale height. 

We adopt the following prescription for the parameter A c 
as a function of stellar mass: 



This is for the stages of core hydrogen burning. During core 
helium burning (CHeB), the value A c = 0.5 is used for all 
stellar masses. 

Overshooting at the lower boundary of convective en- 
velopes is also considered. The value of A c = 0.25 (see Alongi 
et al. 1991, for a description of the formalism) is adopted for 
stars with 0.7 < (M/M Q ) < 2.0. For M > 2.0 M Q a value of 
A c = 0.7 is assumed. 

Mass loss by stellar winds. Perhaps the most uncertain as- 
pect of the evolution of Pop-Ill stars is mass loss by stellar 
winds. Two typical situations can be envisaged: (i) Mass loss 
from massive stars during their whole evolutionary history, and 
(ii) mass loss during the RGB and AGB phases of low and 
intermediate-mass stars. 

Of course, the first question is whether our knowledge for 
normal-metallicity stars can be extrapolated and safely used for 
Pop-Ill stars. This depends on the physical nature of the mech- 
anism driving mass-loss. In massive stars of normal metallicity 
(no matter whether low or high) mass loss at high effective tem- 
perature is due to radiation pressure on resonance lines of ions 
like CIV, NV, SilV, etc., whereas that at low temperatures is 
likely caused by radiation pressure on dust grains and/or due 
to pulsational properties (see Chiosi & Maeder 1986; Chiosi 
et al. 1992; Maeder & Conti 1994; Feast 1991; Kudritzki 1998; 
for recent reviews). Furthermore, the body of observational and 
theoretical information allows us to establish empirical rela- 
tions according to which the mass loss rates are seen to de- 
crease with the metallicity. On this basis one could perhaps 
argue that as long as the surface layers keep the original com- 
position mass loss cannot occur. If and when the surface abun- 
dances are contaminated by mixing of heavy elements, mass 
loss by the above mechanisms could be activated. However, the 
above reasoning would immediately fail if other mechanisms 
(metallicity independent) concur to drive mass loss. 

Even for non-zero metallicities, modelling the RGB and 
AGB phases of low- and intermediate-mass stars is usually 
problematic, as we do not have a definitive theoretical scenario 
to account for mass loss during these stages. Moreover, the 
question arises as to whether the mass-loss prescriptions nor- 
mally used for metal-poor globular cluster giants (e.g. Reimers 
1975), or for pulsating AGB stars (e.g. Vassiliadis & Wood 
1993) can be applied to Pop-Ill stars. Anyway, if any of these 
empirical relations is actually used in zero-metallicity models, 
very little mass is found to be lost during the RGB and AGB 
phases. This justifies to some extent the assumption of evolu- 
tion at constant mass, that we adopt in the present paper. 

However, in Sect. |we explore the possibility that the more 
massive and evolved Pop-Ill stars may become gravitationally 
unbound and hence lose mass. These stars will be the subject 
of a forthcoming paper (Marigo et al. 2001, in preparation). 
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5. General properties of the stellar tracks 
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Table 1. Nuclear lifetimes as a function of the initial stellar 
mass. 
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Fig. 1. Zero-metallicity stellar tracks in the HR diagram, for 
initial masses in the range 0.7 M & - 100 Mq. Evolutionary 
calculations are carried out at constant mass. For the sake of 
simplicity, the He-burning tracks of low-mass (M < 1.1 Mq) 
are not plotted. See also Fig. ||. 

5.1. Mass ranges, evolutionary stages, and lifetimes 

Figure |l] shows the whole set of zero-metallicity tracks at con- 
stant mass in the L — T c g diagram. The initial masses range 
from 0.7 Af Q to 100 M . 

For the sake of clarity in the presentation of the results, 
we adopt the standard convention of sub-grouping the models 
as a function of the initial mass, by defining two mass limits, 
namely AIhcF and M up . 

- Low-mass stars, with M < A/hcF, suffer He-ignition in an 
electron-degenerate core at the tip of the Red Giant Branch 
(RGB). 

- Intermediate-mass stars, with A/hcF < M < M up , ignite 
He in a non-degenerate core, but later develop an electron- 
degenerate C-O core and experience the Asymptotic Giant 
Branch (AGB) phase. 

- Massive stars, with M > M up , are able to ignite central 
carbon in non-degenerate conditions and proceed through 
the sequence of advanced nuclear burnings. 

Our models are evolved from the ZAMS, at constant mass. 
The evolution through the whole H- and He-burning phases is 
followed in detail. The calculations are carried out up to the ini- 
tial stages of the TP- AGB phase in intermediate- and low-mass 
models (a small number of thermal pulses have been followed), 
or the onset of carbon ignition in the helium-exhausted core for 
the most massive models. 

For low-mass stars, the evolution is interrupted at the stage 
of He-flash in the electron-degenerate hydrogen-exhausted 
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core. The evolution is then re-started from a Zero Age Horizon- 
tal Branch (ZAHB) model having the same core mass M^ cF 
and envelope chemical composition as the last RGB model (see 
also Sect. 5.3). 

Nuclear lifetimes of the H- and He-burning phases as a 
function of the stellar mass are reported in Table [j]. 

Let us now discuss the main features of the models here 
presented. It should be recalled that many aspects of the evo- 
lution of zero-metallicity stars are already described by Cassisi 
& Castellani(1993; hereinafter also CC93). In order to avoid 
redundant discussion, the present analysis is more focused on 
particular issues and/or open questions which, to our opinion, 
still deserve further investigation. 
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5.2. The evolution on the H-R diagram 

Let us now discuss the main evolutionary features of zero- 
metallicity stars also in comparison to those of models with 
Z ^ 0. The Z = stellar tracks in the H-R diagram are dis- 
played in Fig. ^, where we mark the onset/end of major core 
burning phases. 

First, we can see that, as expected, the Z = tracks in 
the log L — log T c ff diagram are systematically hotter and more 
luminous. The latter feature results in shorter nuclear lifetimes 
(given in Table [j]) with respect to those of solar-composition 
models. 

Other obvious differences refer to (i) some morphological 
features of the tracks, and (ii) the relative excursion in in effec- 
tive temperature characterising the nuclear H- and He-burning 
phases. 

Concerning point (i), we notice that well-defined loops may 
show up on the Z = tracks (e.g. 1 Mq and 3 Mq models in 
Fig. ||). The physical interpretation of these unusual features 
will be discussed in Sect. ^. 

Concerning point (ii), it is worth making the following re- 
marks. As we can see, the locus of points corresponding to the 
termination of the main sequence (squares) describes a path 
in the logT c ff — logL plane, which significantly bends to- 
wards lower effective temperatures at increasing stellar mass 



(i.e. M <; 20 Mq). This can be explained as due to the gradual 
decrease of the mass exponent a in the mass-luminosity rela- 
tion (i.e. L oc M a ), as pointed out long ago by Stothers (1966). 

Comparing the evolution of 1 Mq models with Z = 
and Z — 0.0004, we see that once they experience powerful 
He-ignition at the tip of the RGB, the subsequent quiescent 
core He-burning occurs either close to the Hayashi line (Z = 
0.0004), or to an HB structure blueward of the RR Lyrae insta- 
bility strip (extending from logT c ff ~ 3.75 to logics ~ 3.85; 
see Bono et al. 1995) when Z = 0. 

Moreover, in stars with M 5M©, the onset of central 
He-burning (circles) occurs soon after hydrogen exhaustion in 
the core, and proceeds as these stars are evolving towards lower 
effective temperatures at nearly constant luminosities. To this 
regard, a feature of particular interest is that the Z = mod- 
els of intermediate mass (1.1 Mq A/ ; < 7 — 8 Mq) start 
and complete their core He-burning phase always in the bluest 
parts of the tracks, i.e. quite far from their Hayashi lines. This 
implies, for instance, that Z = intermediate-mass stars would 
never appear as either red-clump stars or blue-loop stars. 

For even more massive models, the locus of points marking 
the end of the He-burning phase (starred symbols in Fig. |j) 
also presents a systematic bending towards cooler regions at 
increasing stellar masses. In particular, we find that massive 
models with 7OM ^ M Ss 100M Q may reach their Hayashi 
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Fig. 3. Evolution of central conditions 
(i.e. density and temperature) for the 
whole set of Z = models. Approx- 
imate relations are used to determine 
the boundaries (thick straight lines) of 
the regions in which the equation of 
state is expected to be dominated by 
different pressure components, i.e. de- 
generate gas, perfect gas, and radia- 
tion pressure. Labelled lines indicate 
point loci as a function of the stel- 
lar mass corresponding to the onset of 
central H-burning (a); onset of the 3- 
a reaction and hence the CNO-cycle 
(b); end of core H-burning (c); on- 
set of central He-burning (d); end of 
core He-burning (e); energy balance 
(f) between carbon burning and neu- 
trino losses (e c = e^). A few values of 
the initial stellar masses are indicated 
(in M Q ). 
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lines already during the helium-burning phase, ramaining there 
until central carbon ignition. 

In general, it should be remarked that whether or not a star 
spends some part of its evolution on (or in proximity to) the 
Hayashi line may have important implications. In fact, the cor- 
responding development of a deeply extended convective enve- 
lope may cause dredge-up events, with consequent changes in 
the surface chemical composition. This point will be discussed 
in Sect. |} 

5.3. The evolution on the T c — p c diagram 

In addition to the log L — log T c g plane, another powerful diag- 
nostic of the stellar structure is the central temperature - central 
density diagram, as it illustrates the state of the gas in the in- 
nermost regions. Figure [}] displays the evolutionary behaviour 
in this plane by the models here considered. Lines are drawn to 
mark the approximate boundaries of the regions in which the 
equation of state is dominated by perfect gas, degenerate elec- 
tron gas, and radiation pressure. Point loci as a function of the 
stellar mass are shown for relevant evolutionary stages. 

One can notice, in this figure, some typical behaviours that 
depend mainly on the equation of state, and that are found in 
stellar models of any metallicity. For instance, in correspon- 
dence to the stages of gravitational contraction of stellar cores 



between nuclear burnings, the evolution of the central condi- 

1/3 

tions in the domain of ideal gas can be described by T c oc p c , 
as expected from homology considerations (e.g. Kippenhahn & 
Weigert 1990). 

However, the T c — p c diagram of zero-metallicity stars con- 
tain several unusual features worthy of note, such as: the pecu- 
liar loci marking the onset of central H-burning (line a) and the 
onset of the 3-a reaction and the CNO-cycle (line b), as well as 
some short-lived periods of core expansion and/or heating (that 
are evident, for instance, soon above line b), and the unusually 
low mass limit (1.1 Mq) for stars to become degenerate before 
He-ignition. Each one of these features will be discussed below. 

5.4. Up to the onset of the 3-a reaction 

Owing to the complete lack of metals (and hence CNO ele- 
ments) during the "initial" stages of central H-burning, the only 
energy sources able to satisfy the condition of thermal equilib- 
rium are (i) gravitational contraction, and (ii) nuclear burning 
via the p-p chains. In low-mass stars, the p-p chain soon pro- 
vides the nuclear energy necessary to slow down gravitational 
contraction, whereas at increasing mass this energy is reached 
only at higher and higher temperatures. 

The onset of the p-p reactions as a function of mass is very 
inclined in the T c — p c diagram (line a in Fig. ||). Due to the 
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well known rather weak dependence of the p-p reaction rates 
on temperature (i.e. e pp oc T 4 ~ 6 ), the contraction of the central 
regions can easily proceed, so that high temperatures may be 
then attained during the H-burning phase. 

Once the central temperatures have increased up to typi- 
cal values of logT ~ 8, the ignition of the 3-a reaction takes 
place. This occurrence marks a fundamental event, as it is rep- 
resents the first significant production of metals in these stars. 
As a consequence, the synthesis of primary 12 C leads to the ac- 
tivation of the CNO-cycle, which then starts providing nuclear 
energy in competition with the p-p reactions (line b in Fig. ||). 

The onset of the 3-a reaction occurs at earlier and earlier 
stages at increasing stellar mass (see Fig. Q). This occurs es- 
sentially because the 3-a reaction requires much higher tem- 
peratures than the p-p chain. In the lowest mass models, the 
first production of 12 C takes place towards the very end of the 
H-burning phase, and may not even occur for M 0.8 M & . 
In this latter case, all central hydrogen is entirely burnt via the 
p-p chain. At a sufficiently high mass (M > 20M@ in our 
models), the 3-a reaction ignites even before the p-p reactions 
have slowed down the initial stellar contraction. In these mod- 
els then, H-burning simply proceeds via the CNO-cycle, with- 
out any significant phase of central burning via the p-p chain. 

The above features explain the striking change in the slope 
in the curve corresponding to the onset of the H-burning phase 
as a function of the stellar mass (line a in Fig. fy. For mod- 
els with M ^ 20Af Q , H-burning starts above a nearly verti- 
cal line in the T c — p c diagram (i.e. at varying T c and for an 
almost constant p c ), that is characteristic of ZAMS stars that 
burn hydrogen predominantly through the p-p chain. For mod- 
els with M <; 20 Mq, H-burning starts above a line of modest 
slope (i.e. a modest increase of T c corresponds to a substantial 
decrease of p c ), that simply traces the minimum T c and p c nec- 
essary to ignite the 3-a reaction (and hence the CNO-cycle) in 
a core with Y ~ 0.23. This also explains why line (a) simply 
merges with line (b) at M <; 20 M . 

5.5. Soon after the onset of the 3-a reaction 

Subsequent to the first carbon production, as the energy contri- 
bution of the CNO reactions increases (at the expense of the p-p 
reactions) the central regions are forced to expand, with conse- 
quent decrease of the central density (and central temperature 
in the most massive models). This leads to the behaviour of the 
tracks in Fig. ^|, around the line (b). 

The activation of the CNO-cycle may result in the appear- 
ance of some peculiar features in the H-R diagram of low-mass 
models (i.e. loops; see Sect. ||) and, given its quite large tem- 
perature dependence (i.e. £cno °c T 14 ~ 22 ), it also favours the 
development or a larger extension of core convection. For in- 
stance, in the case of the models with M ~ 2 — 5A/ Q displayed 
in Fig. U after the initial recession of the convective core, cen- 
tral H-burning via the p-p reactions occurs in radiative con- 
ditions. Then, as the CNO cycle turns on, central regions be- 
come convectively unstable again. This effect is anticipated at 
increasing stellar mass, and eventually in more massive mod- 
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Fig. 4. The central abundance of hydrogen when the fractional 
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Fig. 5. Evolution of core convection during central H-burning 
for a few selected models with masses (in Mq) as indicated. 

els (M > 5Af Q ) the transition of the core from convective to 
radiative and back to convective does not occur at all. 

5.6. The critical masses: AfucF and M up 

From the location of the non-degenerate/degenerate boundary 
and of the curves for He- and C-ignition shown in Fig. ||, it is 
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Fig. 6. Critical masses A/h f and M up as a function of the 
metallicity according to the predictions of Girardi et al. (2000) 
for values of Z ^ 0, and those presently derived for Z = 0. The 
results of Cassisi & Castellani (1993; and references therein) 
are shown for comparison. 

possible to estimate the critical masses MhcF and M up . Figure 
H illustrates the predicted trend of the mass limits Mh f and 
M up as a function of the metallicity, combining the results by 
Girardi et al. (2000) and those of the present work for Z = 0. 

Concerning Mh 6 f, this is found to decrease with the metal- 
licity, dropping down to ~ 1.1 Mq for Z = 0. This is the com- 
bined result of the larger convective cores during the MS phase 
at decreasing metallicities when the CNO cycle is the dominant 
energy source (up to Z ~ 10~ 4 , see CC93) and, at even lower 
Z, of the higher temperatures reached in the centre when H- 
shell burning mainly occurs via the p-p reactions (as is the case 
for the Z = models). 

Concerning M up , its trend is not monotonic as it first de- 
creases with metallicity, reaches a minimum, and finally in- 
creases again up to ^ 7 Mq - 8 Mq when Z = 0. This be- 
haviour is essentially controlled by the mass of the He-core left 
at the end of the MS. For CNO-dominated H-burning, lower 
metallicities correspond to more concentrated energy sources 
and larger convective cores, which explains the initial decrease 
of M up . As soon as the energy contribution from the p-p re- 
actions becomes competitive with the CNO cycle, the central 
burning regions on the MS are more extended, while the sizes 
of the convective cores tend to be smaller. This explains the 
presence of the minimum and the subsequent increase of M up 
at decreasing Z. 

Both behaviours of Mh f and M up are in agreement with 
the results of CC93. Their larger values of Mh c f and M up at 
given metallicity essentially reflect the different treatment of 



convective boundaries (i.e. Schwarzschild criterion) with re- 
spect to ours (i.e. overshoot scheme; see Sect. Q). 

5.7. Growth of core convection during the He-burning phase 

Another remarkable evolutionary feature is found in stars with 
masses M ~ 1.2 - 2.5 M©, and M £ 30M Q . During the 
central He-burning phase, the outer boundary of the core over- 
shooting zone is located very close to the bottom of the H- 
burning shell. 

Initially, we found that, if the convective core was left to 
grow, it eventually reaches the H-shell and engulfs some H- 
rich material, which is rapidly burnt in the core via the CNO 
cycle. This causes a flash that expands the core, so that cen- 
tral He-burning weakens and the convective core recedes tem- 
porarily (in mass). After the flash has occurred, the convective 
core starts growing again, possibly approaching the bottom of 
the H-shell. This leads to quasi-periodic H-flashes, the related 
growth/recession of the core resembling some sort of breathing 
convection. 

However, this picture is likely to be physically unsound be- 
cause of the the treatment adopted in the evolution code. In fact, 
the H-flash is caused by the fact that the engulfed hydrogen is 
first mixed throughout the convective core, and then burnt ac- 
cording to nuclear reaction rates which are mass-averaged all 
over the convective layers. This procedure is a good approxi- 
mation as long as the convective lifetimes are shorter than the 
nuclear lifetimes, which is usually the case for CNO-burning 
occurring in "ordinary" conditions. 

But at the very high temperatures (~ 10 8 K) reached at the 
bottom the H-burning shell, the nuclear lifetimes involved in 
the CNO cycle (on the order of hours, days) can be comparable 
or even shorter than the typical convective timescales (on the 
order of months). Under these conditions, a correct approach 
would require a time-dependent solution scheme which couples 
simultaneously nucleosynthesis and mixing (see e.g. Schlattl et 
al. 2001) 

Leaving the full analysis of this point to a future investiga- 
tion, we make the reasonable assumption in the present study 
that hydrogen burns locally at the bottom of the H-shell, before 
a complete mixing by core convection can occur. Technically, 
this translates into the condition that the maximum allowed ex- 
tension of core convection is set by the bottom of the H-burning 
shell. 

Figure [7] illustrates the evolution of the convective and 
burning regions of the 2.2 Mq model. We can notice the ap- 
proach of core convection towards the inner location of the H- 
burning shell, during the He-burning phase (third panel from 
left). We can also notice the remarkable thickness (in mass) of 
the H-burning shell, which extends deeply inward during both 
shell H-burning and core He-burning phases. This is due to the 
fact that hydrogen in the shell is mainly burnt via the p-p chain, 
which due to its low temperature sensitivity, allows quite high 
temperatures to be attained at the bottom of the shell. At an age 
of 5.1 x 10 8 yr, our numerical treatment starts to limit the ex- 
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Fig. 7. Kippenhahn diagrams showing the location (in fractional mass) of the burning and convective zones during the evolution 
of the 2.2 M Q , Z = model. 



Fig. 8. Physical properties of the 1 
M© model with Z = at the first 
onset of the 3-a reaction. Top-left 
panel: central abundance of CNO el- 
ements. Bottom-left panel: fractional 
energy contributions from nuclear 
burnings, where L, Lh> ^cno* an d 
L pp refer, respectively, to the surface 
stellar luminosity, the rate of energy 
production via H-burning, CNO cy- 
cle, and p-p chain. Top-right panel: 
evolution of surface quantities as in- 
dicated. Bottom-right panel: evolu- 
tion of central density and tempera- 
ture. 
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tension of the core convection. In this way, the evolution can be 
quietly followed up to the end of He-burning. 



6. Low- and intermediate-mass models 



6.1. The first activation of the CNO cycle 



We have shown in Sect. 5.4 that the onset of the 3-a reaction 
in low-mass models possibly occurs near the end of the central 
H-burning phase. This circumstance is marked by short-lived 
loops in the H-R diagram (see Figs. @ and |). We can better 
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Fig. 9. H-R tracks of Z = low-mass models extending from 
the ZAHB up to the beginning of the TP-AGB phase. 

analyse this feature looking at Fig. g, which displays how rel- 
evant properties of the 1 Mq model vary during the brief time 
interval in which the loop shows up. 

As soon as a CNO abundance of about 10~ 10 '~ 9 is built at 
the centre of the star, a sudden spike arises in the rate of en- 
ergy generation by nuclear H-burning. This occurs when the 
growing efficiency of the CNO cycle exceeds that of the p-p 
reactions (Lcno/L > L pp /L), which have dominated the en- 
ergetics of the star up to then. The central conditions react to 
this sudden increase of energy production: the innermost re- 
gions expand so that both T c and p c decrease. 

6.2. The evolution on the RGB 

As already pointed out by D'Antona (1982) and CC93, a 
prominent feature of zero-metallicity low-mass models is that, 
for given stellar mass, the tip of the RGB is reached with a 
larger core mass and much fainter luminosity with respect to 
models with non-zero metallicity (see, for instance, the 1 Mq 
tracks for Z = and Z = 0.0004 in Fig. g). 

The peculiar behaviour of the luminosity as a function of 
the core mass (M c — L relation) displayed by RGB models 
with Z = can be explained as the combined effect of the 
lower mean molecular weight (/x), and the dominant nuclear en- 
ergy source (i.e. p-p chain). As is shown by Marigo (2000), the 
M c — L relation at Z = can be remarkably well reproduced 
by the homology prediction L oc (jjlM c ) with S ~ 4.5, the 
value of the exponent being essentially determined by the typi- 
cal temperature dependence of the p-p reactions. For dominat- 
ing CNO cycle, we would get 5^7 — 8, which well accounts 
for the luminosity evolution of RGB models with Z > 0. 



Table 2. Characteristics of low-mass stars at the He-flash. 



M 

(Mq) 


M c HcF 

(Mq) 


A 77ZAHB 

^^RGBt 

(10 49 erg) 


A yHcF 


0.7 


0.5015 


2.573 


-0.0325 


0.8 


0.4991 


2.549 


-0.0322 


0.9 


0.4951 


2.497 


-0.0315 


1.0 


0.4882 


2.400 


-0.0303 


1.1 


0.3683 


0.989 


-0.0100 



6.3. The evolution on the HB 

In principle, as a low-mass star settles on the ZAHB, it starts to 
quiescently burn helium in the core with an abundance some- 
what lower than the original one, since a small amount of he- 
lium, Ay c HeF , has been already converted into carbon during 
the flash at the tip of the RGB. 

We have computed this quantity for each evolutionary 
track, by simply assuming that the flash provides just the en- 
ergy necessary to lift the electron degeneracy from the core. In 
practice, this energy can be evaluated as the difference in total 
energy between the last RGB-tip configuration and the initial 
ZAHB one, 

a £rg™ Then ' the 

mass of burnt helium, and 
hence Al^. HcF , follows straightforwardly from the Q-value of 
3-a reactions. Table g presents some relevant characteristics of 
the models at the stages of He-flash and ZAHB. We can notice 
that the predicted fractional abundance of helium is burnt in the 
flash amounts to few percent (typically 3 %). 

The subsequent evolution on the HB (see Fig. |9|) is quite 
similar to that of Pop-II models. However, as already pointed 
out by Cassisi et al. (1996), all low-mass Z — stars burn their 
central He at quite high effective temperatures (with log T e g <; 
3.90), so that a negligible fraction of RR Lyrae pulsators on 
the ZAHB is expected. This aspect is further commented on 
Sect. g. 

6.4. The AGB phase and the occurrence of thermal pulses 

The evolutionary properties of metal deficient AGB stars have 
been analysed in the past by Fujimoto et al. (1984), Chieffi & 
Tornambe (1984), Cassisi et al. (1996), and more recently by 
Dommguez et al. (1999). A crucial question addressed in these 
works is about the possible occurrence of the He-shell flashes 
(thermal pulses), which are usually found in AGB models of 
various (but not zero) metallicities. 

According to the semi-analytical investigation of Fujimoto 
et al. (1984) on the stability of nuclear burning shells, ther- 
mal pulses are conditional on the existence of thresholds in 
core mass M c , and CNO abundance .^cno m the envelope. 
To summarise, He-shell flashes are expected to occur when i) 
M c < 0.73 M Q for any value of Z CNO , and ii) M c > 0.73 Mq 
if .Zcno - > ^CNO> tn e critical lower limit depending on M c . 
The dependence on Zcno can be qualitatively understood con- 
sidering that in the absence (or strong deficit) of CNO elements 
in the envelope, the H-burning shell itself produces new carbon 
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Fig. 10. Evolutionary properties dur- 
ing the first thermal instabilities of 
the He-burning shell, for two Z = 
AGB models with initial masses of 
2.5 M , and 5.0 M s . See text for fur- 
ther explanation. 



via the 3-a reaction, so that the CNO cycle can operate. It fol- 
lows that the H- and He-burning shells are not energetically de- 
coupled, and proceed outward in mass at the same rate (see also 
Chieffi & Tornambe 1984). Since He-rich material is not accu- 
mulated in the inert intershell buffer, but rather steadily burnt 
by the underlying He-shell, the occurrence of thermal pulses is 
thus prevented. 

In this study we present the results for two intermediate- 
mass models, with initial zero metallicity and masses of 
2.5 M Q and 5 M Q , which are made evolve through the initial 
stages of the AGB phase, including a few episodes of thermal 
instabilities of the He-burning shell. Whether such events may 
be ascribed to true He-shell flashes is discussed below. 

As the 2.5 Mq and 5 M Q models reach their Hayashi tracks 
(E-AGB phase) after central helium exhaustion, envelope con- 
vection progressively moves inward in both cases, and extends 
into regions which have experienced earlier nucleosynthesis. 
The so-called "second dredge-up" takes place, but with a sub- 
stantial difference between the two models. This can be appre- 
ciated by looking at Table Bl which presents the predicted sur- 
face chemical compositions after the second dredge-up event 
for all our models. 



In the 2.5 M Q model the envelope penetrates into the chem- 
ical profile left by the recession of the convective core during 
the MS, and into regions previously affected by the p-p chain, 
so that the surface enrichment is essentially in helium with no 
trace of heavier metals (see also Fig. [7]for a similar case). In the 
5 Mq model the envelope penetration proceeds further, reach- 
ing down to the former location of the H-burning shell (which 
temporarily extinguishes), so that newly synthesized CNO ele- 
ments are also dredged-up to the surface. After the occurrence 
of the second dredge-up, the evolution of these AGB models 
proceeds quite differently. 

As we can see from Fig. [HJ, the 2.5 M model experiences 
rather weak fluctuations in the surface luminosity, reflecting 
the alternation between the H- and He-burning shells in pro- 
viding nuclear energy to the star. The long-term behaviour of 
the surface properties indicates the model is evolving down on 
its Hayashi track, i.e. to decreasing luminosity and increasing 
effective temperature. As already mentioned, the CNO surface 
abundance is zero, so that the H-burning shell must itself pro- 
duce the necessary 12 C to sustain the CNO-cycle. The latter 
is actually the dominant energy source during most part of the 
AGB evolution followed by our calculations (bottom-left panel 
of Fig. p. 

The 5 M Q model suffers well-defined He-shell flashes, 
which drive the features usually found in the light curves of 



M/Mq 


H 


3 He 




4 He 


12 C 




13 C 




14 N 




15 N 




16q 




17 o 




18q 




22 Ne 




24+25 Mg 


20 Ne 




Initial: 


















































all 


0.770 


Z.40 1U 


'S 


U.Z3U 


U 




u 




U 




u 




U 




f\ 

u 




u 




() 









n 
U 




After the first dredj 


;e-up: 














































0.7 


0.769 


8.63 10" 


-5 


0.231 




















































0.8 


0.768 


1.4010" 


_4 


0.232 




















































0.9 


0.768 


1.5410" 


-4 


a oqo 
U.Z32 


U 




A 
U 




A 

u 




A 

u 




A 

\) 




A 

U 




A 

U 




A 
\) 









A 
U 




1.0 


0.768 


1.1710" 


-4 




U 




a 
U 




U 




A 

u 




A 

U 




A 
\) 




A 

u 




A 

u 









A 

u 




1.1 


0.769 


8.29 10" 


-5 


n n 1 


A 

U 




A 

u 




u 




A 
\) 




A 
\) 




A 

u 




A 

u 




A 

u 









A 
U 




After the second dredge-up: 














































1.2 


0.769 


6.49 10" 


-4 


0.231 




















































1.3 


0.767 


4.89 10" 


-4 


0.233 




















































1.4 


0.769 


4.7210" 


-4 


0.231 




















































1.5 


0.767 


3.45 10" 


-4 


0.233 




















































1.6 


0.767 


3.1610" 


-4 


0.232 




















































1.7 


0.767 


2.75 10" 


-4 


0.233 




















































1.8 


0.769 


2.64 10" 


-4 


0.231 




















































1.9 


0.768 


2.1910" 


-4 


0.232 




















































2.0 


0.768 


2.02 10" 


-4 


0.232 




















































2.1 


0.704 


1.0710" 


-4 


0.296 




















































2.2 


0.697 


9.5910" 


-5 


0.303 




















































2.3 


0.691 


8.1310" 


-5 


0.309 




















































2.4 


0.690 


6.8210" 


-5 


0.310 




















































2.5 


0.681 


5.9710" 


-5 


0.319 




















































2.7 


0.669 


5.0410" 


-5 


0.331 


4.1710" 


-21 


1.3410" 


21 


4.95 10" 


19 


2.02 10" 


-23 


1.28 10" 


20 


1.3310" 


■22 


7.1610" 


-26 


2.43 10" 


ss 


1.6210" 


-50 


7.70 10" 


ss 


3.0 


0.660 


4.1010" 


-5 


0.340 


1.13 10" 


-18 


3.6210" 


19 


1.2810" 


Ifi 


5.1910" 


21 


3.15 10" 


-18 


2.82 10" 


■20 


1.6010" 


-23 


4.25 10" 


-32 


6.91 10" 


-46 


3.1910" 




3.5 


0.649 


2.99 10" 


-5 


0.351 


8.06 10" 


-16 


2.5910" 


16 


8.5810" 


14 


3.4610" 


-18 


1.97 10" 


■IS 


1.41 10" 


17 


8.69 10" 


-21 


2.06 10" 


-28 


1.8610" 




6.1910" 


■30 


4.0 


0.644 


2.29 10" 


-5 


0.356 


6.33 10" 


-14 


2.03 10" 


14 


6.1710" 


12 


2.46 10" 


-16 


1.31 10" 


■13 


7.41 10" 


-16 


5.23 10" 


-19 


5.33 10" 


-26 


7.4210" 


-37 


4.13 10" 


-27 


5.0 


0.637 


1.56 10" 


-5 


0.363 


3.5810 


-9 


7.81 10" 


12 


6.06 10" 


■10 


1.57 10" 


-14 


4.54 10" 


12 


3.9010" 


14 


5.11 10" 


-15 


5.22 10" 


-19 


6.68 10" 


-28 


1.1410" 


23 


6.0 


0.632 


1.11 10" 


-5 


0.368 


2.63 10 


-7 


1.1910" 


11 


9.35 10" 


■10 


2.41 10" 


-14 


7.3610" 


11 


5.91 10" 


14 


4.67 10" 


-13 


2.96 10" 


Ifi 


1.6610" 


- 23 


4.73 10" 


21 


6.5 


0.626 


9.25 10" 


-6 


0.374 


7.07 10 


-7 


1.4610" 


11 


1.1310- 


-9 


2.90 10" 


-14 


2.99 10" 


■10 


7.37 10" 


14 


2.21 10" 


-12 


1.95 10" 


■15 


2.63 10" 


-22 


4.09 10" 


■20 


7.0 


0.668 


8.74 10" 


-6 


0.332 


1.26 10" 


-13 


4.04 10" 


14 


1.0910" 


11 


4.3010" 


16 


2.08 10" 


■13 


8.57 10" 


Ifi 


7.51 10" 


-19 


6.02 10" 


26 


2.75 10" 


-36 


9.03 10" 


■27 


8.0 


0.681 


7.49 10" 


-6 


0.319 


9.66 10" 


-14 


5.68 10" 


17 


8.08 10" 


-12 


6.3210" 


-19 


1.47 10" 


-13 


1.41 10" 


-18 


1.15 10" 


-21 


9.09 10" 


■30 


6.05 10" 


-42 


4.43 10" 


■31 


8.3 


0.687 


7.0010" 


-6 


0.313 


1.1910" 


-14 


3.81 10" 


15 


1.01 10" 


12 


3.98 10" 


-17 


1.88 10" 


14 


7.3910" 


-17 


6.71 10" 


-20 


1.95 10" 


■27 


2.1210" 


-38 


2.3610" 


■28 


9.0 


0.770 


4.3010" 


-5 


0.230 




















































9.5 


0.770 


4.3010" 


-5 


0.230 




















































10.0 


0.770 


4.25 10" 


-5 


0.230 




















































12.0 


0.770 


3.63 10" 


-5 


0.230 




















































15.0 


0.770 


2.47 10" 


-5 


0.230 




















































20.0 


0.770 


1.6610" 


-5 


0.230 




















































30.0 


0.770 


1.0210" 


.5 


0.230 




















































50.0 


0.770 


5.95 10" 


-6 


0.230 




















































70.0 


0.713 


1.0610" 


-6 


0.287 


7.65 10" 


-12 


2.45 10" 


12 


6.39 10" 


■10 


2.53 10" 


-14 


1.3010" 


11 


5.29 10" 


14 


5.3410" 


-17 


5.11 10" 


24 


2.08 10" 


-14 


2.58 10" 


14 


100.0 


0.640 


6.39 ID- 


-7 


0.360 


2.27 10" 


-11 


7.25 10" 


12 


1.87 10 


-9 


7.21 10" 


-14 


3.53 10" 


-11 


1.33 10" 


■13 


1.33 10" 


-16 


2.03 lo- 


23 


1.4210" 


-32 


1.13 lo- 


23 



14 



P. Marigo et al.: Zero-metallicity stars 



AGB stars. Namely, over a pulse cycle, the luminosity under- 
goes the initial post-flash rapid dip and peak, which is followed 
by the low-luminosity dip, and the final recovering of the lu- 
minosity level that characterizes the quiescent shell H-burning 
(see, for instance, Boothroyd & Sackmann 1988). In this case, 
the long-term evolution is characterised by the increase of the 
surface luminosity and slight decrease of the effective tempera- 
ture, i.e. the model is climbing on the AGB. The envelope CNO 
abundance (by mass) left after the completion of the second 
dredge-up is Z C no = 4.198 x 1(T 9 . 

As for the 2.5 M® model, the CNO cycle provides most 
of the stellar energy during the stages of quiescent H-burning 
(bottom-right panel of Fig. |l0|). However, it should be noticed 
that a fraction of the CNO nuclei involved as catalysts in the 
CNO cycle are not brought up to the surface by the second 
dredge-up, but are synthesized in situ by the H-burning shell 
via the 3 — a process. This can be seen by comparing the total 
CNO abundance (by number) in the envelope, Zqno, an d that, 
for instance, at the point of maximum nuclear efficiency within 
the H-burning shell, ^'no- ^ there were not production of 12 C 
inside the shell, we should have Zqno = ^cno> whereas we 
find Zcno < ^cno- with ^cNO ranging from - 2 x 1(T 8 
to ~ 6.5 x 10~ 8 at subsequent stages immediately preceding a 
thermal pulse. 

Let us now check whether our results can be recovered 
within the Fujimoto scheme. The 2.5 Mq and 5 Mq models 
have core masses of about 0.69 Mq and 0.87 Mq, and Zcno 
in the envelope equal to and ~ 4.198 x 10 -9 , respectively. 
According to Fujimoto et al. (1984) the lowest mass model is 
expected to suffer thermal pulses, since M c < M™* (regard- 
less of 2cno). whereas the most massive one should not, since 
it is located inside the region of the plane M c — Zcno where 
thermal pulses are prohibited (see their figure 5). Thus, our re- 
sults would seem to disagree with the Fujimoto scheme. This 
is perhaps not surprising as the predictions by Fujimoto et al. 
(1984) are derived from a semi-analytical method, in which the 
input physics - adopted for the integrations of deep envelope 
structures - may differ from those presently used in our stellar 
code, as well as for other model prescriptions, such as the ini- 
tial helium abundance (Y = 0.25 in Fujimoto et al., Y = 0.23 
in our models). Other differences may also arise from the use 
of analytical approximations in the Fujimoto procedure. 

However, even if our results do not strictly obey the quan- 
titative thresholds in M c and Zqno as pointed out by Fujimoto 
et al. (1984), a qualitative agreement in the behaviour of AGB 
models with initial zero metallicity can be found. In particular, 
we confirm the prediction that the occurrence of thermal pulses 
in massive AGB models is conditional on the pre-existence of 
some CNO abundance in the envelope, that in the case of our 
5Mq model is supplied by the second dredge-up during the 
E-AGB. This finding is also supported by the evolutionary cal- 
culations of zero-metallicity intermediate-mass stars by Chieffi 
& Tornambe (1984) and Dommguez et al. (1999). 

The strength of thermal instabilities should also depend on 
the 2cno va hie in the envelope, as this quantity controls the 
degree of thermal coupling between the H- and He-burning 



Table 4. Core masses and central C/O ratio at the stage of C- 
ignition. 
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shell. In models with initial zero metallicity, we expect that the 
more Zqno is dredged-up, the less 12 C is synthesised by the 
3-a reaction in H-shell during a interpulse period, hence the 
stronger the He-shell flash should be. 

7. Massive models 

7.1. Core masses 

In Figure [H] and Table Q we report the masses of the He and 
C-O cores of the zero-metallicity models with M > M up ~ 
7 Mq, which are expected to end their evolution as supernovae. 

These quantities represent fundamental information for 
both hydrodynamic calculations of the explosion event (see, 
for instance, Umeda et al. 2000), and synthetic derivation of 
the supernova yields (see the procedure described by Portinari 
et al. 1998). 

Comparing the results for Z = with those for other metal- 
licities, we can notice that: i) for initial masses Ss 25 — 30 Mq 
the core masses are almost independent of metallicity; ii) the 
major differences show up at higher masses, being largely de- 
termined by the effect of mass loss that is expected to be less 
efficient at lower Z (and not applied to the Z = models); iii) 
in the case of the Z = models evolved at constant mass, both 
A/hc and Mqo follow a nearly linear relation with the stellar 
mass. 

7.2. The Eddington critical luminosity 

We also check whether the most massive models may become 
gravitationally unbound, i.e. the rate of energy outflow at the 
surface exceeds the corresponding Eddington luminosity 



L K = 



inGcM 



(5) 



where M is the stellar mass, and k the opacity of the surface 
layer. 
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Fig. 11. Masses of the He core (Mh c ) and C-O core (Mco) at the central carbon ignition in zero-metallicity massive models, 
compared with predictions for Z = 0.02, and Z = 10 -4 . 



In Fig. [l2| we compare the evolutionary paths in the H-R 
diagram for selected models of different initial masses with the 
loci obtained from Eq. (Q), i.e. calculating the Eddington lumi- 
nosity of the photosphere for each value of the effective tem- 
perature along the evolutionary sequence. 

From this simple test, it turns out that, in contrast with the 
50 A/© model, the 70 Mq and 100 M© models may achieve 
super-Eddington luminosities (typically at logT e g ~ 3.9) on 
their way to the Hayashi line, towards the end of the He- 
burning phase. Correspondingly, the effective acceleration]] at 
the surface becomes negative, due to the increase of the surface 
radiative opacity (see Fig.|l3|). 

From this point onward the predicted evolution should not 
be considered reliable, as the adopted assumption of hydro- 
static equilibrium does not hold any longer for the outermost 
layers. Likely, these stars would start losing mass from their 
surface. 

For the sake of simplicity, we do not attempt to include any 
prescription for mass-loss driven by super-Eddington luminosi- 
ties in the most massive stars in this work. This, and other stel- 
lar winds driving mechanisms, will be analysed in more detail 
in a future paper dedicated to the evolution of zero-metallicity 
stars with mass loss (Marigo et al. 2001, in preparation). 




Fig. 12. Evolutionary tracks in the H-R diagram (thick lines) 
for selected massive models with initial masses (in A/©) as in- 
dicated. Thin-line curves labelled with Le show the behaviour 
of the Eddington luminosity calculated at the photospheric 
layer with Eq. (|5|). 

8. Surface chemical changes 

We will now discuss the changes in the surface chemical com- 
position due to convective dredge-up episodes. Usually these 

1 The effective acceleration is defined as GM/R 2 — (cr/c) k.T^ s , 
i.e. the difference between the gravitational and radiative accelera- 



16 



P. Marigo et al.: Zero-metallicity stars 




Fig. 13. Effective acceleration at the photosphere, defined as 
the difference between the inward gravitational acceleration, 
ffgrav, and the outward radiative acceleration, g ra d. A sequence 
of models are shown up to the point in which the net acceler- 
ation becomes negative, as the surface luminosity exceeds the 
Eddington limit. Correspondingly, the evolution of the photo- 
spheric opacity, k, is shown. 



result from the penetration of the convective envelope into stel- 
lar layers which have previously been modified by nucleosyn- 
thetic processes. In this respect it is relevant to recall that the 
deepest inward extension of envelope convection occurs as a 
star reaches close to its Hayashi line. Keeping this in mind and 
looking Table [3] we can make the following considerations. 

Between the end of the H-burning phase and the central 
He-ignition, the only models that evolve to the red part of the 
H-R diagram and settle onto the giant branches (closely ap- 
proximating the location of the Hayashi tracks) are those of 
low-mass (0.7 M© <, M & 1.1 Mq). At this stage envelope 
convection starts to proceed inward, but it does not extend sig- 
nificanlty into the chemical profile left by the radiative burn- 
ing core. As can be noticed in Table ||, this first dredge-up 
episode in low-mass stars changes the surface chemical abun- 
dance of helium by negligible amounts. For more massive mod- 
els (M > 1.2 Mq) the onset of the He-burning phase takes 
place far away from the Hayashi tracks. As a consequence, the 
first notable feature of zero-metallicity models is that the first 
dredge-up essentially does not take place. 

As far as the second dredge-up is concerned, we can notice 
from Fig. ^| that after central helium exhaustion the evolution 
for all models with M < 70 M© is characterised by a red- 
ward excursion in the H-R diagram. However, the minimum 
value of the effective temperature attained - either at the start 
of the AGB phase or at central carbon ignition - has a non- 
monotonic dependence on the stellar mass. Specifically, stars 
with initial masses 0.7 M© ^ M ^ 8 A/© are able to ap- 
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Fig. 14. Surface abundances of 4 He and CNO elements as a 
function of the stellar mass, for the last computed model. 



proach their Hayashi tracks, then becoming red giants. In this 
mass range the second dredge-up is found to actually occur 
in models with 2.1 Mq M S$ 8.3 Mq, the envelope pen- 
etration being insufficient at lower masses. Moreover, we can 
notice from Table | that for 2.1 M© <S M < 2.1 Mq only 
helium increases its abundance in the envelope, whereas for 
2.7 Mq Si M Ss 8.3 Mq some primary CNO is also dredged- 
up to the surface. 

It is worth remarking that further changes in the sur- 
face chemical abundances may be caused by the occurrence 
of the third dredge-up at thermal pulses (if they take place, 



see Sect. 6.4) during the TP- AGB phase (see, for instance, 



Dommguez et al. 1999). Actually, in our calculations we do not 
find any evidence of this process, but this is more likely due to 
the small number of He-shell flashes that have been followed 
for each evolutionary sequence. 
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More massive stars (8.3 M Q < M < 50 M ) met the nec- 
essary conditions for carbon ignition in the core at higher effec- 
tive temperatures, i.e. before entering the giant branch phase. 
As expected, these models do not suffer any dredge-up episode, 
i.e. their surface layers preserve the original composition. 

Finally, stars with 70 M© < M < 100 M Q settle on their 
Hayashi tracks already during the core He-burning phase, and 
there they remain until the onset of carbon burning. The sec- 
ond dredge-up takes place enriching the surface composition 
both in helium and CNO elements (see Table ||). As already 



discussed in Sect. [7.2 
nous models 



the results for these massive and lumi- 
made to evolve at constant mass - should be 
taken with some caution. Our analysis has in fact revealed that 
these models may become dynamically unstable against radi- 
ation pressure towards the end of core helium burning. From 
that point, stellar winds driven by super-Eddington luminosi- 
ties will indirectly affect the surface chemical composition by 
i) reducing the stellar mass, and ii) influencing the efficiency 
of the second dredge-up. To this regard, a more careful study 
is going to be carried out in a forthcoming work (Marigo et al. 
2001, in preparation). 

In conclusion, owing to the particular evolutionary features 
of the zero-metallicity stars, we find that significant changes 
in the chemical composition of the envelope - during the evo- 
lutionary phases here considered - may occur within two well 
defined mass ranges, i.e. 2.1 M < M < 8 M Q , and 70 M Q < 
M < 100 M Q . Our results are fully illustrated in Fig. [l4|. 



OB 

o 




4.4 



Fig. 15. Instability strips in the log L — log T c g diagram. Theo- 
retical edges for (Z = 0, Y = 0.23) are drawn assuming pulsa- 
tion occurs in the fundamental mode (solid line), first overtone 
(dashed line), and second overtone (dotted line). The shaded 
area shows the expected location of RR Lyrae according to 
Bono et al. (1995). See text for more details. 



9. Might we ever detect primordial stars ? 

Zero-metallicity stars with initial masses M < 0.8 M , if ever 
born, should be still alive, as their nuclear lifetimes are longer 
than the age of the universe (see Table [[]). However, to date the 
hypothesis of the possible formation of primordial low-mass 
stars has not been proven observationally. In order to devise a 
test for the existence of such stars, one should single out ob- 
servable distinctive features which can point unambiguously 
to their existence. For instance, the absence of metal lines in 
the stellar spectra would immediately give us the confirmation, 
but this circumstance is hard to meet. In fact, the interaction 
of long-lived first stars with a chemically enriched interstellar 
medium would have probably polluted their primordial surface 
compositions with newly synthesized metals. 

Another peculiar feature of low-mass zero-metallicity stars 
(with M ~ 1 Mq) is the development of a loops in the H-R 
diagram, corresponding to the first activation of the CNO cycle 
as soon as some primary 



12 C is produced by the 3-a reaction 



(see Sect. 5.4). Unfortunately, the probability of detecting such 



loop in a primordial simple stellar population is negligible, as 
this feature is quite short-lived (a few 10 6 yr), hence scarcely 
populated. 

Finally, distinctive observational properties of low-mass 
zero-metallicity stars may be derived from their loci in the HR- 
diagram: they are slightly hotter than the Pop-II stars of even 
the lowest metallicity. For instance, Cassisi et al. (1996) find 
that most of the evolution of Z = 10~ 10 HB stars occurs to the 



left (higher T e g) of the RR Lyrae instability strip, thus imply- 
ing a very low probability of finding RR Lyrae pulsators among 
an hypothetical population of these stars. 

However, we notice that high effective temperatures are 
also met at the turn-offs of low-mass stars. This might open 
the possibility that, at ages > 10 Gyr, the main-sequence of 
Z = populations falls inside the faintest part of the instability 
strip which contains (in order of decreasing luminosity) three 
classes of variables, namely: Cepheids, RR Lyrae and S Scuti 
stars. In such a case, the possibility arises that a large number of 
short-period, low-amplitude pulsators might be present among 
Pop-Hi stars. 

To investigate this possibility, the results of a preliminary 
analysis of pulsation and evolution of Z = stars are shown in 
Figs. [l5] and[l(| We have made use of the large data set of the- 
oretical models of Cepheid variables calculated by Chiosi et al. 
(1993), corresponding to luminosities \ogL/L & > 2, various 
values of the stellar mass (3M @ <, M <, 12M Q ), and dif- 
ferent metal and helium abundances (Z = 0.02, 0.008, 0.004; 
Y =0.23, 0.24). This data set has been extended towards 
lower luminosities by means of a preliminary stability anal- 
ysis, specifically performed for zero-metallicity models with 
M = IMq. The linear nonadiabatic pulsation calculations 
have been carried out with the aid of the same code as in Chiosi 
et al. (1993; and references therein), but incorporating the low 
temperature opacities of Alexander and Ferguson (1994). The 
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new calculations cover the region of the HR diagram delimited 
by 0.3 < log(£/L ) < 1.3 and 3.90 > logT cfr ^ 3.73. 

Combining the earlier pulsation models with the new ones, 
we are now able to determine the edges of the instability strip 
from logL/Z/Q ~ 4.5 down to log L/L & ~ 0, as illustrated in 
Fig. [l5| Analogously to Chiosi et al. (1993), we have also de- 
rived several parameterised fitting formulas (given in Table ||) 
as a function of helium Y and metal content Z. They describe, 
for different pulsation modes, the location of the instability 
strips and other relevant relations between period, luminosity, 
colour, and mass. The instability strips displayed in Fig. |l5| are 
obtained from the aforementioned analytical relations for the 
red and blue edge setting Z = and Y = 0.23. 

The main interesting points are the following. At luminosi- 
ties log L/Lq > 2 the instability strip crosses the evolutionary 
tracks of primordial stars which should already be dead and 
hence undetectable. Towards lower luminosities the theoretical 
location of the instability strip is consistent with the predictions 
for RR Lyrae variables according to Bono et al. (1995) for fixed 
helium and metal content (Y = 0.24, Z = 0.001). The shaded 
region corresponds to the maximum extension of the strip (to 
the red and to the blue), combining the results of Bono et al. 
(1995) for both fundamental and first overtone pulsators. We 
can notice that the RR Lyrae instability strip is almost twice as 
extended to the red as our instability strip at the same luminosi- 
ties. It is worth remarking here that, whereas the blue edge is 
well determined by our stability analysis, the red edge is not 
(for more details see Chiosi et al. 1993). As a conservative ap- 
proach, we define the red edge to coincide with the maximum 
of the growth rate. In this way, we are certainly underestimating 
the width of the instability strip. 

Anyhow, at the typical luminosities of RR Lyrae, the insta- 
bility strip for the Z = models would cover the red portions 
of the HBs described by the zero-metallicity low-mass tracks 
plotted in Fig. |l5](see Cassisi et al. 1996, for a similar result). 
In other words, pulsation should set in during a period of very 
rapid evolution, i.e. only after these stars have burnt most of 
their central helium in the bluer regions close to the ZAHB. 
This again would make the observability almost out of reach. 

Finally, going to even fainter luminosities (logi/L© Ss 1) 
it turns out that zero-metallicity models with mass as low as 
0.9 Mq are expected to exhibit pulsational properties during 
the main-sequence phase. At the luminosities and effective 
temperatures under consideration (see Fig. |l5|), the pulsation 
periods can range from about 0.6 to 0.03 days. The reliability 
of these provisional results is supported to some extent by the 
fact that the location on the H-R diagram and the period deter- 
minations are consistent with the observational data of 5 Scuti 
stars - a rich class of faint variables with amplitudes ranging 
from hardly detectable to several tenths of a magnitude - which 
are usually present among the Blue Stragglers of many Galactic 
globular clusters. A sample of observed 5 Scuti stars (Nemec 
& Mateo 1990; McNamara 2000) is shown in the M v versus 
B — V diagram of Figure [lf^, and compared to the theoret- 
ical instability strip. It is interesting to notice that, applying 
the fundamental-mode period-mass-luminosity-colour relation 



-4 



-2 




Fig. 16. Instability strips in the My — (B — V) diagram, with 
the same notation as in Fig. |l5[ Isochrones of ages 15.8 Gyr and 
10.0 Gyr are shown for Z = (continuous line) and Z = 10~ 4 
(from Girardi et al. 1996; dot-dashed line). Points correspond 
to observed S Scuti stars in various galactic globular clusters 
and in the Carina DSph galaxy (Nemec & Mateo 1990, filled 
cirles; McNamara 2000, starred symbols). 

presented in Table ||, to the S Scuti data (Nemec & Mateo 1990; 
McNamara 2000) and assuming Z = 0.001, we can immedi- 
ately estimate the pulsational masses for these variables, with a 
typical range between 1.2 and 1.8 Mq. 

A direct comparison with observations is presented in 
Fig. [lf^, showing isochrones representative of Pop-II (Z = 
10 -4 ) and Pop-Hi (Z = 0; see the Appendix) stellar popu- 
lations, for two values of the age, i.e. 10.0 Gry and 15.8 Gyr. 
These latter define the best range for the age of Galactic glob- 
ular clusters, as recently reviewed by Carretta et al. (2000). 

We can see that in none of the four cases displayed, the 
turn-off region is located inside the instability strip, except for 
the Z = 0, 10.0 Gyr old isochrone, which marginally inter- 
sects the red edge of the strip. However, we note that, since our 
calculations probably underestimate the red extension of the 
instability strip (see earlier in this section), there is still some 
possibility that the Z = isochrones for ages within the range 
under consideration will enter significantly into the red part of 
the instability strip. If this is indeed the case, then it implies that 
a large number of primordial pulsators could be present in the 
present era. Such stars would presumably be similar to 5 Scuti 
variables but with lower masses at a given colour (since Z = 
stars near the main-sequence are bluer and more luminous than 
solar metallicity stars of the same mass). They would also have 
longer periods at a given luminosity and colour (since pulsation 
period varies approximately as M~ - 5 ). 
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Table 5. Fitting relations for pulsational models. Coefficients are given for the fundamental mode, first and second overtone. 
Respective periods (P , Pi, and P 2 ) are expressed in days, masses and luminosities are in solar units, helium (Y) and metal (Z) 
abundances are by mass. Standard notation of the quantities is used. 



Mode 


A 


B C D E F 


G 


H I 




Blue Edge: lo 


g T SE =A + BxY + CxZ + Dx (logM) + E x (logL) + F x (logL) 2 


' +GxY(logL) 


2 + HxZ(logL) 2 


Po 
Pi 
P 2 


3.928 
3.962 
3.936 


2.985E-2 -5.021E-1 2.986E-2 -5.130E-2 -1.311E-3 
-2.0117E-1 -6.903E-1 8.774E-2 -1.298E-2 -1.945E-2 
4.816E-2 -8.542E-1 8.414E-2 -3.018E-2 -1.515E-2 


5.915E-3 
3.073E-2 
1.758E-2 


-2.319E-2 
-2.658E-2 
6.488E-2 




Red Edge: logT fl , B =A + BxY + CxZ + Dx (logM) + E x (logL) + F x (logL) 2 


+ GxY(logL) 2 


+ HxZ(logL) 2 


Po 
Pi 
P 2 


3.879 
3.901 
3.876 


5.948E-2 2.274E-1 6.708E-2 -4.646E-2 -3.762E-3 
-7.233E-3 5.884E-2 6.731E-2 -4.071E-2 -6.632E-3 
1.134E-1 -5.327E-1 2.224E-2 -3.824E-2 -3.246E-3 


-4.628E-4 
6.754E-3 
2.557E-3 


-1.511E-1 
-1.177E-1 
3.747E-2 


Blue Edge: (B-V)ss 


= A + BxY + CxZ + Dx(Mv) + Ex(YxMv) +Fx(ZxM v ) + Gx(M 2 ^) + Hx(YxM 2 / ) + Ix(ZxM 2 / ) 


Po 
Pi 
P 2 


3.102E-1 
-2.215E-1 
1.062E-1 


-1.925E-1 5.349 -5.217E-2 4.720E-2 -2.590E-1 

1.653E 8.027 -3.594E-1 1.406 2.767 
3.175E-1 4.482 -2.750E-1 1.004 2.314 


6.541E-3 
-2.355E-2 
-2.932E-2 


-1.019E-2 1.516E-1 
1.524E-1 7.409E-1 
1.664E-1 7.820E-1 


Blue edge: (V-I) BB : 


= A + BxY + CxZ + Dx(M v ) + Ex(YxMy) + Fx(ZxM ¥ ) + Gx(M^) + Hx(YxM 2 v ) + Ix(ZxM^) 


Po 
Pi 
P 2 


5.029E-1 
8.104E-2 
3.184E-1 


-4.736E-1 2.331 -3.438E-2 -9.280E-2 -1.857E-1 
8.586E-1 3.731 -2.736E-1 9.622E-1 1.408 
-1.931E-1 1.809 -2.387E-1 7.655E-1 2.243 


6.263E-3 
-1.644E-2 
-2.424E-2 


-2.05 1E-2 3.047E-2 
1.071E-1 3.687E-1 
1.377E-1 6.635E-1 


Red edge: (B-V)i?s 


= A + Bx Y + CxZ + Dx(Mv) + Ex(Y*M v ) + Fx(ZxM F ) + Gx(M 2 / ) + Hx(YxM^) + Ix(ZxM^) 


Po 
Pi 
P 2 


3.349E-1 
1.341E-1 
3.258E-1 


1.407E-1 4.187 -1.405E-1 3.711E-1 -1.638 
8.584E-1 6.649 -2.564E-1 8.899E-1 1.708 
7.263E-2 4.307 -1.777E-1 5.641E-1 2.764E-1 


-5.235E-3 
-2.054E-2 
-1.799E-2 


5.933E-2 3.714E-1 
1.237E-1 9.079E-1 
9.780E-2 2.353E-1 


Red edge: (Y-l) RE ■ 


= A + BxY + CxZ + Dx(Mi/) + Ex(Y*Mi/) + Fx(ZxM v ) + Gx(M 2 v ) + Hx(YxMy) + Ix(ZxM 2 v ) 


Po 
Pi 
P 2 


5.525E-1 
3.827E-1 
5.315E-1 


-1.107E-1 3.013E-1 -9.441E-2 1.655E-1 -1.415 
4.601E-1 2.334 -1.970E-1 6.056E-1 1.088 
-1.852E-1 6.837E-1 -1.495E-1 4.114E-1 1.909E-2 


-1.833E-3 
-1.605E-2 
-1.790E-2 


3.032E-2 1.165E-1 
8.835E-2 5.258E-1 
8.645E-2 2.669E-2 


logP = A + B x (logL) + C x (logL) 2 + D x logM + E x logT oft 


Po 
Pi 
P 2 


12.150 
10.256 
10.002 


0.659 3.000E-2 -0.723 -3.567 
0.683 2.269E-2 -0.671 -3.112 
0.692 1.978E-2 -0.685 -3.070 






lo 


gP = AxlogM + B + CxZ + DxMy + Ex(ZxMy) + Fx(B-V) + Gx(Z 2 x(B-V)) + Hx(B-V) 2 


' +Ix(Z 2 x(B-V) 2 ) 


Po 

Pi 
p 2 


-6.656E-1 
-6.072E-1 
-5.760E-1 


-3.771E-1 -7.788 -3.411E-1 4.206E-1 6.834E-1 
-5.719E-1 -8.367 -3.254E-1 2.656E-1 8.998E-1 
-7.067E-1 -1.016E+1 -3.099E-1 -4.138E-1 1.035 


5.034E+2 
4.089E+2 
1.092E+3 


1.095E-1 -2.927E+2 
-7.956E-2 -1.571E+2 
-1.451E-1 -1.345E+3 




logP = A x logM + B+CxZ + DxMv+Ex(ZxMv) + Fx (V-I) + G x (Z 2 x (V-I)) + H x (V-I) 2 


+ Ix(Z 2 x(V-I) 2 ) 


Po 

Pi 
p 2 


-6.864E-1 
-6.103E-1 
-5.661E-1 


-3.636E-1 -2.742 -3.347E-1 -8.353E-1 2.243E-1 
-5.764E-1 -5.814 -3.159E-1 -1.017 4.879E-1 
-7.197E-1 -6.262 -3.007E-1 -1.382 6.468E-1 


1.317E+2 
2.659E+2 
6.420E+2 


5.397E-1 -1.859E+2 
3.382E-1 -2.217E+2 
1.995E-1 -8.657E+2 



Might we ever detect them ? 
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Appendix A: Tables of evolutionary tracks 

The data tables for the present evolutionary tracks are avail- 
able only in electronic format, either upon request to the au- 
thors, or by accessing a WWW site containing a complete 
da ta-base (including addit ional data and the future extensions) 
at frittp://pleiadi. pd.astro.it . 

For each evolutionary track, the corresponding data file is 
organised into El] columns with the following information: 



1. age/yr: stellar age inyr; 

2. logL: logarithm of surface luminosity (in solar units), 

log(L/L Q ); 

3. logTef: logarithm of effective temperature (in K), 

logT cff ; 

4. grav: logarithm of surface gravity ( in cgs units); 

5. logTc: logarithm of central temperature (in K); 

6. logrho: logarithm of central density (in cgs units); 
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7. Xc, Yc: mass fraction of either hydrogen (up to the central 
H-exhaustion) or helium (later stages) in the stellar centre; 

8. Xc_C: mass fraction of carbon in the stellar centre; 

9. Xc_0: mass fraction of oxygen in the stellar centre; 

10. Q_conv: fractional mass of the convective core; 

1 1 . Q_di s c: fractional mass of the first mesh point where the 
chemical composition differs from the surface value; 

12. L_H/L: the luminosity provided by H-burning reactions as 
a fraction of the surface luminosity; 

13. Q1_H: fractional mass of the inner border of the H-rich re- 
gion; 

14. Q2_H: fractional mass of the outer border of the H-burning 
region; 

15. L_He/L: the luminosity provided by He-burning reactions 
as a fraction of the surface luminosity; 

16. Ql_He: fractional mass of the inner border of the He- 
burning region; 

17. Q2_He: fractional mass of the outer border of the He- 
burning region; 

18. L_C/L: the luminosity provided by C-burning reactions as 
a fraction of the surface luminosity; 

19. L_nu/L: the luminosity lost by neutrinos (hence negative) 
as a fraction of the surface luminosity; 

20. Q_Tmax: fractional mass of the point with the highest tem- 
perature inside the star; 

21. stage: label indicating particular evolutionary stages. 

A number of evolutionary stages are indicated along the 
tracks (column namely: the initial evolutionary stages 
(ZAMS or ZAHB), local maxima and minima of L and T e g 
(Te-M, Te-m, L-M, and L-m), the exhaustion of central hy- 
drogen (Xc=0) and helium (Yc=0), the lowest L and highest 
T e g during the He-burning of intermediate-mass stars (Bhe 
and LpT, respectively), the base and tip of the first ascent of the 
red giant branch (Brg and Tip, respectively), the maximum L 
immediately preceding a thermal pulse (ltp), and the onset 
of C-burning (Cb). These stages delimit characteristic evolu- 
tionary phases, and can be useful for the derivation of physical 
quantities (as e.g. typical lifetimes) as a function of either mass 
or metallicity. Notice that some of these evolutionary stages 
may be absent from particular tracks, depending on the precise 
value of stellar mass and metallicity. 



Appendix B: Tables of isochrones 

From the tracks presented in this paper, we have constructed 
isochrones adopting the same algorithm of "equivalent evolu- 
tionary points" as used in Bertelli et al. (1994) and Girardi et al. 
(2000). The initial point of each isochrone is a 0.7 M Q model in 
the main sequence. The terminal stage of the isochrones corre- 
sponds either to stars at the beginning of the TP-AGB for ages 
<; 10 8 yr, or to stars at C-ignition for ^ 10 8 yr. The youngest 
isochrones (up to 6 Myr) are basically limited to the main se- 
quence stages, since at these ages the inclusion of more evolved 
stars would require tracks for M 100 A/©, that will be pre- 
sented in a forthcoming paper. 
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Fig. B.l. Theoretical isochrones in the HR diagram for the 
initial composition [Z = 0,Y = 0.23]. The age range goes 
from log(t/yr) = 6.3 to 10.2, at equally spaced intervals of 
A log t = 0.3. In all isochrones, the main sequence is complete 
down to 0.7 M Q . 



Complete tables with the isochrones can be obtained 
upon request to the authors, or through the WWW site 
frittp://pleiadi. pd.astro.it They cover the complete age range 
from about 2 Myr to 16 Gyr (6.3 < log(*/yr) < 10.25). 
Isochrones are provided at A log t = 0.05 intervals; this means 
that any two consecutive isochrones differ by only 12 percent 
in their ages. 



Figure B.l shows some of these isochrones on the HRD, 
sampled at age intervals of A log t = 0.3. 

Theoretical luminosities and effective temperatures along 
the isochrones are translated to magnitudes and colors by us- 
ing extensive tabulations of bolometric corrections and colors, 
as in Bertelli et al. (1994). The tabulations were obtained from 
convolving the spectral energy distributions contained in the 
[M/H] = -2.0 library of stellar spectra of Kurucz (1992) 
with the response function of several broad-band filters. The re- 
sponse functions are from Buser & Kurucz (1978) for the UBV 
pass-bands, from Bessell (1990) for the R and / Cousins, and 
finally from Bessell & Brett (1988) for the JHK ones. 

For each isochrone table, the layout is as follows: A header 
presents the basic information about the age and metallicity of 
each isochrone. Column 1 presents the logarithm of the age in 
yr; columns 2 and 3 the initial and actual stellar masses, re- 
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spectively. We recall that the initial mass is the useful quantity 
for population synthesis calculations, since together with the 
initial mass function it determines the relative number of stars 
in different sections of the isochrones. Then follow the loga- 
rithms of surface luminosity (column 4), effective temperature 
(column 5), and surface gravity (column 6). From columns 7 to 
15, we have the sequence of absolute magnitudes, starting with 
the bolometric one and following those in the U BVRIJHK 
pass-bands. In the last column (16), the indefinite integral over 
the initial mass M of the initial mass function (IMF) by num- 
ber, i.e. 



is presented, for the case of the Salpeter IMF, (f>(M) = AM~ a , 
with a = 2.35. When we assume a normalization constant of 
A = 1, FLUM is simply given by FLUM= A'/ 1_Q /(1 - a). This 
is a useful quantity since the difference between any two values 
of FLUM is proportional to the number of stars located in the 
corresponding mass interval. It is worth remarking that FLUM 
relations can be derived for alternative choices of the IMF, by 
using the values of the initial mass we present in the column 2 
of the isochrone tables. 

We also provide summary tables containing basic informa- 
tion for the most significant stages along the isochrones. The 
evolutionary stages listed are, in sequence: 

- TO: the turn-off point, i.e. the point of highest T g during 
the core-H burning phase. 

- If present, Te-m and Te-M signal the coldest and hottest 
points, respectively, of stars in the overall contraction phase 
at the end of core— H burning. 

- RGBb and RGBt: the base and tip of the RGB, respectively. 

- BHeb: the beginning of the CHeB phase. 

- If present, Te-m and Te-M signal the coldest and hottest 
points, respectively, for CHeB stars. For the youngest 
isochrones, Te-M represents the maximum extension of the 
Cepheid loop. 

- EHeb: the end of the CHeB phase. 

- In the oldest isochrones, L-M and L-m limit the luminosity 
range of early-AGB stars; this interval corresponds to the 
clump of early-AGB stars in colour-magnitude diagrams. 

- 1TP: the beginning of the thermally pulsing AGB phase. 

- Cb: the stage of C-ignition in the cases the AGB phase do 
not occur. 

In addition, we provide tables with the integrated broad- 
band colours of single-burst stellar populations. Such tables are 
derived by integrating the stellar luminosities, weighted by the 
predicted number of stars in each bin, along the isochrones. 




(B.l) 



